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Abstract. The curvature tensor and the scalar curvature are computed in the space of positive 
definite real matrices endowed by the Kubo-Mori inner product as a Riemannian metric. 
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1. Introduction 

The state space of a finite quantum system is identified with the set of positive scmidefinite 
matrices of trace 1. The set of all strictly positive definite matrices of trace 1 becomes naturally 
a differentiable manifold and the Kubo-Mori scalar product defines a Riemannian structure on 
it. Reference [4] tells about the relation of this metric to the von Neumann entropy functional. 

The objective of the paper is to compute the scalar curvature in the Riemannian geometry of 
the Kubo-Mori scalar product. Actually, we consider the space of real density matrices which is a 
geodetic submamifold in the space of complex density matrices. Our study is strongly motivated 
by the conjectures formulated in [3] and [4]. It was conjectured that the scalar curvature takes its 
maximum when all eigenvalues of the density matrix are equal, and more generally the scalar cur- 
vature is monotone with respect to the majorization relation of matrices. Although we obtain an 
cxplicitc formula for the scalar curvature, the conjecture remains unproven. Nevertheless, a huge 
number of numerical examples are still supporting the conjecture. The method of computation 
of this paper is inspired by [2]. 

When this paper was nearly finished we recieved the preprint [1] where the scalar curvature is 
computed for the Kubo-Mori metric in the complex case by a different method. 

2. The Kubo-Mori metric on the space of positive definite matrices 

2.1. The setup. Let S = S(n) be the space of all real selfadjoint (n x n)-matriccs, 6> + = S+(n) 
be the open subspace of positive definite matrices. Then S+ is a manifold with tangent bundle 
TS+ = S+ x S. We shall consider the following Riemannian metric on S+, where D £ S+ and 
X,Y e T D S+ = S: 

POO 

G D {X,Y) = / Tr ({D + t^XiD + ty^dt. 
Jo 
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Because 

Tr ((£> + t^XiD + t^Y) | < nr 2 \\X\\ \\Y\\ 

the integral is finite. We shall identify S with its dual S* by the standard (i.e., Hilbert-Schmidt) 
inner product (X, Y) = Tr (XY). Then we can view the Riemannian metric G also as a mapping 



_D • J £><->+ — <-> — > J £><->+ — <-> = <->> 

poo 

G D (X)= / (D + ^-^^ + t)" 1 ^, 
Jo 

which is symmetric with respect to (A", Y) = Tr XY. (Note that Gd is the Frechet derivative 
of log.D.) Now let D e <S + and choose a basis of R™ such that D = ^ Aj£jj is diagonal, where 
(Eij) is the usual system of matrix units, then the selfadjoint matrices 

F kl = E kl + E lk (k < I) 

are a complete system of eigenvectors of Gd ■ S — > S. This means that (F i} )i<i<j< n is an 
orthogonal basis of (TpS, Gd) with 

for (»,.?) ,6 (M) 

G D (Fij,F k i) = { 2m, 3 for i = fc < j = / 
4TOii for z = j = k = I 

/°°(a, + 1) - \x t + tr dt = =: Wfei . 

The expression m k i is a symmetric function of the eigenvalues X k and A;. In fact l/m k i is the 
logarithmic mean of X k and A; . This implies that 

m kl = -!- whenever A fe = A;, 

in particular, m kk — l/X k . Note that GoiFij) — rriijFij for all i < j. 
Another symmetric expression 



where 



/ (Ai + t)- 1 (X j + t)-\X k + t)- 1 dt = TO f ~ ^ =: m yfe 
Jo *i ~ A & 



will appear below. The identity 



1 / m kki m k g \ _ ^ 
m k i \m kk m u 



is easily computed and will be used later. 

2.2. The Christoffel symbol. Since we have a global chart we can express the Levi-Civita 
connection by one Christoffel symbol: 

(V iV )\ D = d V (D)4(D) - r D (£(D), r,(D)) 

where £, r\ : S+ —* S are smooth vector fields. The Christoffel symbol is then given by 

G D (T D (X,Y),Z) = \dG(D)(Z)(X,Y) - \dG(D)(X)(Z,Y) - \dG{D){Y){X,Z), 

where the derivative of the metric 

dG(D)(Z)(X, Y) = J Tr {-{D + t)~ 1 Z{D + t)~ 1 X{D + t)^Y- 

-(D + ty 1 x{D + ty x z{D + dt 
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is visibly symmetric in the entries X, Y. 
The Christoffel form is given by 



G D (T D (X,Y)) = -\dG(D){X){Y) 

= \ J (D + ty 1 (X{D + t^Y + Y(D + t)- 1 *) (D + t)- 1 dt, 
T D (X, ) = -\G- D \dG{D){X) 



Since 



G^{X) = f D U XD 1 - U du, 
Jo 



we can express the Christoffel form as an integral formula. The derivative is 

dT(D)(X)(Y, Z) = -\d{G- l ){D){X).dG{D){Y){Z) + ±G^ d 2 G(D)(X,Y)(Z). 

2.3. When D = J2i \Eu is diagonal, then 

dG(D){Z){X) = — '^^m i j k (E ii Z Ejj X E kk + EnXEjjZE kk ). 

ijk 

In particular, 

dG(D)(Fij)(F k i) = -SjkmujFu - 6jim ikj F lk - 8 l im i , ki Fj k - S lk mjuFji, 

and 

r D (F 21 ,F kl ) = -5 jk ^F u - 5^F lk - Su^F jk - 6 ik ^F jt 
m a m lk m jk rriji 

2.4. The curvature. The Riemannian curvature R(£,ri)( = (V^V,, — V^V^ — V[£ i7? ])C is then 
determined in terms of the Christoffel form by 

R D (X, Y)Z = -dT{D){X){Y, Z) + dT(D)(Y)(X, Z) + 
+ T D (X, T D {Y,Z))- Y D (Y, Y d (X,Z)). 

If we insert the expressions from 2.2 we get after some computation 

R D (X,Y)Z = \d(G- l )(D){X).dG{D){Y)(Z) - |d(G- 1 )(i?)(y).dG(^)(X)(Z) 
= -\G D \dG{D){X).G D \dG{D){Y){Z) 
+ \G D \dG{D){Y).G D \dG{D){X){Z) 

The Ricci curvature is then given by the following trace 

Ric D (X, Z) = Tr s (Y^R D (X,Y)Z), 

and the scalar curvature is 

Scal(L>) = Tr $ (X i — > G D 1 .Ric I) (X, )). 

Next we compute the traces in a concrete basis. Let A s be an orthonormal basis with respect 
to the inner product (X, Y) = Tr (XY) on S. Then 

Kic D (X, Z) = Y,( r d(X, A s )Z, A s ) 

s 

and 

Scal(D) = ^(G^ 1 . Bic D (A t , ), A t ) = ]T (Ric D (A 4 , ),G D 1 A t ) 
t t 

= Mc D (A t , A*) = ^2^2(R D (A t ,A a )Go 1 A t ,A a ) . 
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3. The submanifold of normalized matrices 

3.1. The submanifold of trace 1 matrices. We consider the affine submanifold of 5+ of all 
positive definite real selfadjoint matrices with trace 1 and its tangent bundle: 

S 1 = {DeS+:Tr(D) = l}, 

TSi = Si x S , where S = {X E S : Tr (£>) = 0}. 

Lemma. 

(1) The unit normal field n along the submanifold Si with respect to the Riemannian metric 
G from 2.1 is given by n(D) = (D,D). 

(2) The G-orthonormal projection Pd ■ S = TdS + — > To Si = So is given by Pd{X) = 
X - Tr (X)n{D) = X - Tr (X)D for D E Si and X e S). 

(3) The Christoffel form for the covariant derivative V 1 of the induced Riemannian metric 
Si, G 1 is given by 

T^iXX) = P D T D (X,Y) = Y D (X,Y) - Tr (T D (X,Y )).£>, D E Si, X,Y E So, 
and i/ie second fundamental form is given by 

S D : TdSx x T^Si = S x S -» R 

&(x,r) = Tr(r D (x,r)) 

/■ oc 

= / Tr({D + t)- 1 X{D + t)- 1 Y 
Jo 

- \D{D + t)- 2 X{D + t)- l Y 

- \D{D + t)- 2 Y{D + dt. 

Proof. If X E S commutes with D we get 

G D {X,Y) = Tr ((D + t)- 1 X(D + t)- 1 Y)dt= Tr ((Z) + t)~ 2 XY)dt = 
Jo Jo 

= [- Tr ((£> + t^XY)]^ = Tr (D^XY) 

Thus for F e S we have G D (D, Y) = Tr (Y) = 0. Moreover for D E Si we have G D (D, D) = 
Tr(D) = 1, so (1) follows. The remaining assertions are standard facts from Riemannian geom- 
etry. 

For the explicit expression of the second fundamental form we preceed as follows. For D E Si 
the Weingarten mapping is given by 

Ld -TdSi = So — > So = TdSi, 

L D (X) : = V (AX) n = dn{D).X - T D (X, «(£>)) =X- T D (X, D), 

and the second fundamental form is then given by 

S D :T D Si x TbSi = So x S -» R 

S D (X,Y) : = G D (L D (X),y) =G D (X-r D (X, J D),F) 

= G D {X,Y) + \dG{D){X){D,Y) 

= / Tr ((D + ^-^D + ^Y 
Jo 

- \D{D + t)- 2 X(D + t^Y - \D{D + t)- 2 Y(D + dt. 

Another formula for the second fundamental form is 

(4) S D (X,Y) : = Tr(r d (X,Y))=Tr(-±G D 1 dG(X)(Y)) 

= \\ / Tr (D u (D + t)- 1 X{D + t)- 1 Y{D + t)- 1 D 1 - u + 

Jo Jo 

+ d u (d + t) _1 Y(D + ty 1 x{D + ty 1 D 1 - u ^ dtdu. a 
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3.2. The curvature via the Gaufi equation. The Gaufi equation expresses the curvature R 
of S + and the curvature R 1 of Si for D e Si and X, Y,Z,U£ So by 

G D {R{X,Y)Z, U) = G D (R\X, Y)Z, U) + S D (X, Z)S D (Y, U) - S D (Y, Z)S D (X, U). 

The Ricci curvature of the submanifold S\ is then given by the following trace 

Ric 1 ,^, Z) = Tr 5o (Y -> R D (X,Y)Z) , 

and the scalar curvature is 

Scal 1 ^) =Tr So (X^G^ 1 .Rici,(X, )). 

Next we compute the traces in a concrete basis in case of a diagonal D = J2i KEu G Si. Let 
A s be an orthonormal basis with respect to the inner product Go on So- Then 

Ric 1 , (X, Z) = Y,G D (Rd (X, A S )Z, A s ) 

8 

and 

Scal 1 ^) = Y, Gd{G d \ Ric}, (A t , ), A t ) = £<Rici,(A t) ),A t ) 
t t 

= Y Rici (A t ,A t ) = Y52GD(R 1 D (A t ,A a )A t ,A a ) 

t t s 

= Y(G D {R D {A u A s )A t ,A s ) - S D (A t ,A t )S D (A s ,A s ) - S D {A s ,A t )S D {A t ,A s )) 

4. Computation of the scalar curvature 

Our aim is to have an explicit formula for the scalar curvature Scal 1 (D) in terms of eigenvalues 
the of D, which we assume to be a diagonal matrix. As in the previous section, let A s be an 
orthonormal basis with respect to the inner product Gd on So- We assume that some of the 
basis elements are diagonal (like D) and the others are normalized symmetrized matrix units. 

4.1. The first term. We decompose the sum 

Y j G d (Rd{A u A s )A u A s ) 

t,s 

into three subsums and we compute them separately. First we consider the case when both A t 
and A s are offdiagonal, that is, they arc in the form / ^/2my . 
Offdiagonal-offdiagonal. 

E 4m,-m H Gl ' (flp(fl] ' ft ' )fl] ' Ft ' ) = E -^-{R D (F l0 ,F kl )F l3 ,F kl ) 

+ Y,^^ G D 1 AG)(D).{F kl ).G-jAG(D)(F l3 ){F l0 ),F kl ) 1 

where summation is over i < j and k < I. We continue with the first term and calculate in an 
elementary way: 

- E (Go 1 AG) (D) (Fy ).G~ D \dG{D) (F kl ) {F l3 ),F kl ) 
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For the second term we use dG(D)(Fij)(Fij) — —mujFa — rriijjFjj and get 
+ E Y^( G D-d(G)(D).(F kl ).G^.dG(D)(F ij )(F ij ),F kl ) 

i<j,k<l lJ 

1 x ^ Wluuv^T'uuw . 1 \ A '^'vvw'^uvv . 1 \ A 
9 — ' rn 771 ?tj 9. — ' m m m 9. ^- — ' 



UWW ' ' ' C ■!■ (!■ 



TfX ixdTTI uw Tfl uu 2 717 vw 717 uv 717 vv 2 717 uw 717 vw 717 h 

u<v<w u<v<w u<v<w 



4 — ' m l .m :; A ^ — ' 



Offdiagonal-diagonal. Next we compute the sum 

Y,Gd(Rd{A u A s )A u A s ) 

t,s 

when A t = J2i a \Eu are diagonal, Gd -orthogonal to D, and orthonormalized, and where the A s 
are still offdiagonal. This means that 

, a t a t' 

G D {D,A t ) = ^m ii \ i a\ = ^ j a\ = Q and G D (A t , A t >) = E fn u a\a\ = ^ -i-i- = S t ,f. 

i i ii 

We also have 

dG(D)(F kl )(A t ) = dG(D)(A t )(F kl ) = -(m kkl al + m kll a\)F kl , dG(D)(A t )(A t ,) = -D~ 2 A t A t , , 
dG{D)(F k i)(F k i) = -5 k i(m uk + m kk i)F k i - m kk iF kk - m uk F u , and G D x A t = DA t 

since D, A t , and A t > commute. We get 

E 7^—G D (R D (A t ,F kl )A t ,F kl ) 
t% 2mkl 



Eg^W) ) -dG(D) (F kl )(A t ), F kl ) 



t,k<l 



t,fc<; 



t,k<l 



+ E e^T - (G^ 1 -dG(D) (A t ) (A t ), dG(D)(F k i ) (F k i)) = 
= - ^ ( (-m kk ia k - m k ua\)F kU (-m kk ia k - m kU a\)F kl ) 

t.k<l 

+ ^2 E T~( a p) 2 ^PP' {- m kkiF k k - mmFii)) = 

t,fe</ 8TOfez P A p 

= - E T-2-( m kkia t k +m kll a t l ) 2 
t.k<i TOfei 

+ E 4^ (-«^(4) 2 +rn Hk ^(a\f) =: Q 
t,fe<; » / 

Denoting by Q this seemingly basis dependent quantity we transform the sums in Q as follows: 

E(-H(E(-)-E(->V 

t,fe</ yt,*;,; t,fe=; y 
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Summing for k = I indexes we obtain: 

- E 7^r( m kki4 + m kll a\f = - ^(4) 2 = - 7 E 72 (4) 2 , 

t.k=l t.k 4 t ,k A * 

£ ik (-4«> 2 + "•4<°'> 2 ) - 5 ii^r «> 2 - 1 5 i«> 2 

The two terms turned out to be equal, so 



-2-(m fc fcio fc +m W io,) +2^7— H— +^—( a i) 



We start to deal with the first sum. Let 



V Afc 



Then A, & 1 , b 2 , . . . , b^~^ is an orthonormal basis in R™. We define a linear mapping K. from R™ 
to the space of all real n x n matrices (endowed by the standard Hilbert-Schmidt inner product). 

JCc = y^ — !— [rrikki + m kii y/hci] E k i (c e R"). 

^ ' 777.1.; L J 



Then 



ll^ &t || 2 = E ^2" [ m kkial + m k uaj] 

which is a term in Q. Hence 



t,fc,i m kl 



Since 



and 



||iCA|| 2 = J2 [mkklXk \ mkllXl] = £1 = „■ 



2 



2 



Tr jc*jc = v imi 2 = y 2^k^ + 2 y 

Z 2,/ u I 11 



we have 



E i( m «a4 + m *« a ') 2 = 2 £ + ^ £ x; 1 

t,kd m kl k ,l m kl Z i 



n 2 . 



The other terms of Q are similarly computed as traces. 

E l / mkki , t s 2 , m "fc / t\2\ \ ~> ™fcfci + mm _ \ - m kk i\ k + m k uXi 
TTJ...I \ X,, ^ fc ' \, l ' I ~ 2-~t rn,., ^ 



l hi m kl \ \ k k} A, y u J fa m kl ' fa " m kl 



E 



mkkl + m ku 2 
n 



k,i mkl 
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Finally, we obtain a basis independent expression for Q: 

n _ \ " m kki^k _ _J_ ,_i \ - mkki + rn ku 

Q ~h ^ ~ 16 v * tr *»« 

Diagonal-ofFdiagonal. This case is completely similar (in fact, symmetric) and yields the 
same Q. 

Diagonal-diagonal. Now we compute the sum 

Y,Gd{Rd(A u A s )A u A s ) 

t,s 

when A t and A s are both diagonal, G_D-orthogonal to D, and orthonormalized. We get 
Y j G d {Rd{A u A s )A u A s ) 

t,s 

= -Y,\^ G ^ D )^)- G D 1 - d G{D){A s ){A t \A s ) 

t,s 

+ E \ M G ) ( D ) ( A ' )- G ~d - dG ( D ) ( A t )(At),A s ) = 

t,s 

= -Y,l( G D-dG(D)(A s )(A t ),d(G)(D)(A t )(A s )) 

t,s 

+ J2l( G D 1 -dG(D)(A t )(A t ),d(G)(D)(A s )(A s )) = 

t,s 

= - E \((- D ~ lA t A *> -D- 2 A t A s ) - {-D~\A t f, -D- 2 (A S ) 2 )) = 

t,s 

4.2. The second term. We start the computation of 

-J2S D (A t ,A t )S D (A s ,A s ). 

t.s 

We use first the formula from 3.1.(4), and use also 2.3. 

S D (F t3 ,F kl ) = Tr (-^G-J dG(F t] ){F k i)) 

= Tr (—^^(—SjkmujFu — SjiniikjFik — SumjkiFjk — SikmjuFji)) 

l rr ( e m Uj , c m ikj ^ , r mjki ^ , £ rriju \ 
= \ Tr S jk -Fa + ^7 J -F lk + 8ii^—F jk + kk^—F fl 

V m U m ik m jk rUjl J 

= OjfeOji h djidik h Oiidjfe 1- djfedji 

my m ifc mjk rtiji 

We observe that for (i < j) 7^ (k < I) we get So(Fij, F k i) = 0. Furthermore, for i < j we have 



'33 
TTliki 



ma m 
Sd(Fu, F k k) — 4<5ifc 

o r TP TP \ XX m Ui 1 r r m ifci j r "lit! r r 

oD{fii,-Pkl) = dikdu h dj/djfc h di'Oifc h OjaAz 

ran m ik rriik m a 

= if k < I. 
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First we take summation when both A t and A s are offdiagonal: 
-J2S D (A t ,A t )S D (A s ,A s ) = 

El ( Tnuj | rriijj \ ( rrikki , ram 
dm ■ ■ iji} 

I VJ l- n2 ( n -^ 



i<j,k<l J N JJ ' ^ 



4.^, 16 



Next we assume that A t is diagonal and A s is offdiagonal. 



2 

p 



and 



,j 1'"/ / V m kk m u J 16 

Note that this contribution is symmetric in t and s and we should take it twice. 
Now it is easily follows the value of the sum when both A t and A s are diagonal: 

- ]T S D (A t ,A t )S D (A s ,A s ) = - £ ± . I = -±(„ - I) 2 

t,s t,s 

4.3. The third term. Now we have to deal with 

- J2 S D (A u A s )S D (A s ,A t ) = - VJ S D (A u A s f 

By the formulas from 4.2 we have that it suffices to sum when both A t and A s are diagonal and 
both of them are offdiagonal. Hence 

-Y,Sn{A u A s )S D{ A s ^ 
Tl J j£ x Am kl V ™kk m u J 4 8 



4.4. The scalar curvature formula. 



3 m uvw ^ 1 Tfl vvw Tfl vvu 

A Z — ✓ T 77 m 7TJ 9 ^ — ' 



4 m JUV Tfl vw Tfl'i i } U 2 Tfl vu Tfl vv Tfl 

u<v<w u<v<w 



2 m, wu Tn wv Tn ww 2 m JUU Tft uv Tn uw 

u<v<w u<v<w 

8 ^ ^-mu 8 ^ m ii m jj 

E mjkl^k _ 1 \p , _i \ - m fcfc ; + mfcH 

n 2 (n-l) 2 n(n-l) 2 1 ,, 2 n-1 n 2 (n-l) 2 

(n - l) 2 — 

16 8 4 V ; 4 16 
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Some further simplification: 



8 ^ ro?,m« 8 ^ mf-m^ 16 4^ m?, " 2 ^ A, 



1 m fc fc ; A fc _ _Iy 
2 tT 4 Tf m ^ 



Hence 
R = 



3 v-^ m 



2 

uvw 



4 Tnuv^lvw^lwu 
u<v<w 



j 1 V > TlXyyyjTTiyyxi ^ 1 V > 'H^WWuTTl'WWV ^ 1 V ^ 7Tl uuv 7Tl 
9 ' ^ -m . rrj _ _.rn 9 ' ^ rrj w. .rrj 9 ^ — ^ 



li li 1' ii u if 



2 TTlyu'fTlvv'fW'VW 2 771 wu 77~l'wu 711 ww 2 WluuWluv'ITluw 

u<v<w u<v<w u<v<w 

1 V- m tki X k + m 2 kll Xi 13 ^ 1 s^ mm+mjM n(n - 1) ^ 2 

i6 tr ro « 2 r A * tr 4m « 4 
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